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Abstract. Firstly, we give a short review about the hydrodynamic model and its 
application to the elliptic flow phenomena in relativistic heavy ion collisions. Secondly, 
we show the first approach to construct a unified model for the description of the 
dynamics in relativistic heavy ion collisions. 



1. A Short Review on Elliptic Flow from Hydrodynamic models 

First data reported by the STAR Collaboration at RHIC pQ has a significant meaning 
that the observed large magnitude of elliptic flow for charged hadrons is consistent with 
hydrodynamic predictions [2 . This suggests that large pressure possibly in the partonic 
phase is built at the early stage (r ~ 0.6 fm/c) in Au+Au collisions at ^snn — 130 
and 200 GeV. This situation at RHIC is in contrast to that at lower energies such as 
AGS or SPS where hydrodynamics always overpredicts the data [Hj. Moreover, this also 
suggests that the effects of the viscosity in the QGP phase is remarkably small and 
that the QGP is almost a perfect fluid jl]. Hadronic transport models are very good 
to describe experimental data at lower energies, while they fail to reproduce such large 
values of elliptic flow parameter at RHIC (see, e.g., Ref. [3]). So the importance of 
hydrodynamics is rising in heavy ion physics. After the first STAR data were published 
PP, other groups at RHIC have also obtained the data concerning with flow phenomena 
jSj- To understand these experimental data, hydrodynamic analyses are also performed 
extensively (ZlEj- In this short review, we highlight several results on elliptic flow from 
hydrodynamic calculations. 



1.1. Why elliptic flow? 

Elliptic flow is very sensitive to the degree of secondary interactions |J]. The indicator 
of momentum anisotropy is the second harmonic coefficient of azimuthal distributions 
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Figure 1. Comparison of the hydrodynamic result with the results obtained by a 
partonic transport model. Figure taken from Ref. |14|. 
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For recent progress of higher harmonics, see Refs. ^IJ ^] . The spatial anisotropy 

= (y 2 -x 2 ) s f---e(T,x)d 3 x 

(x 2 + y 2 ) s ' K '"' s fe(r,x)d 3 x 1 ' 

at the initial time is a seed of the momentum anisotropy. Within the hydrodynamic 
picture, pressure gradient along the x-axis is larger than the y-axis in non-central 
collisions. Here the direction of x-axis is parallel to the impact parameter vector. This 
is the driving force for generating t> 2 - The response of the system can be defined by v-ile. 
As will be discussed later, the hydrodynamic response is almost constant v 2 /e ~ 0.2-0.24 
at the SPS and RHIC energies [T3] . 

In ideal hydrodynamics, the mean free path among particles is assumed to be zero. 
So the hydrodynamic prediction of v 2 should be maximum among transport theories. 
We naively expect that the ideal hydrodynamics corresponds to the infinite limit of 
the cross sections among particles in collision terms of the Boltzmann equation. In 
Ref. Molnar and Huovinen compare ^(pt) from hydrodynamics with the one 

from Boltzmann equation. Special attention is paid that the same initial conditions 
and thermodynamic properties are used in comparison. Figure [0 shows the transverse 
momentum dependence of v 2 for a hydrodynamic model and a parton cascade model. 
t>2 (j>r) from Boltzmann simulations is approaching to the hydrodynamic result with 
increasing the cross section among partons. However, even for large cross section 
(47 mb), the result from the Boltzmann simulation is still ~ 30% smaller than the 
hydrodynamic result. This clearly shows the ideal hydrodynamics describes the strongly 
interacting matter. 
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1.2. Basics of Ideal Hydrodynamics 

Here we concentrate our discussions on the ideal hydrodynamic model. For those who 
have interest in the viscous hydrodynamics, see, e.g., Ref. [TH] . 

Hydrodynamic equations are nothing but the energy-momentum conservations 
d p T^ ,u = 0. In the ideal hydrodynamics, the energy- momentum tensor becomes 
T^ u = (e + P)u^u v — Pg^, where e is the energy density, P is the pressure, and 
is the four fluid velocity. When there are conserved quantities such as the baryon 
number or the number of chemically frozen hadrons, one needs to solve the continuity 
equations <9 M nf = together with the hydrodynamic equations. In order to close the 
system of partial differential equations, the equation of state (EoS) P(e, ni) is needed. 
It is the EoS that governs the dynamics of the system in the ideal hydrodynamics. The 
naive applicability conditions of ideal hydrodynamics are that the mean free path among 
the particles is much smaller than the typical size of the system and that the system 
keeps local thermal equilibrium during expansion. From these conditions, one cannot 
use hydrodynamics for initial collisions, final free streaming and high pt 2 GeV/c) 
particles. The quantitative criteria are more complicated in general since these come 
from dynamical aspects of the space-time evolution. 

One usually assumes the initial hydrodynamic fields at initial (or equilibrated) time 
around 1 fm/c. One needs an interface between the pre-thermalisation stage and the 
hydrodynamic stage at the initial time. On the other hand, the system eventually 
breaks up and the assumption of the thermalisation is no longer valid at the later stage. 
This means a prescription of freezeout is needed in the hydrodynamic model when one 
evaluates the particle spectra to be observed by the detectors. Then one needs another 
interface between the hydrodynamic stage and the free streaming stage. In what follows, 
we discuss particularly equation of state, initial condition and freezeout prescription used 
in the literature. 

1.2.1. Equation of State The main ingredient of the hydrodynamic model is the 
equation of state (EoS) for thermalised matter produced in heavy ion collisions. Ideally, 
one uses the EoS taken from the first principle calculations of QCD, namely, lattice QCD 
simulations ^B|- More practically, one can use the resonance gas model for the hadron 
phase and the massless free parton gas for the QGP phase. By matching these two EoS's 
at the critical temperature, one obtains the first order phase transition model with a 
latent heat ~1 GeV/fm 3 . In the mixed phase in this model, the sound velocity c 2 s = ^ 
is vanishing. Recent lattice QCD simulations tell us that the phase transition seems 
to be crossover in vanishing baryonic chemical potential. Although the discontinuity of 
the thermodynamic variables do not exist in the crossover phase transition, it should be 
emphasized that the energy density and the entropy density increase more rapidly than 
the pressure in the vicinity of the phase transition region AT ~ 0.1T C . This also leads 
to very small sound velocity near the phase transition region. Therefore it is very hard 
to find flow observables which distinguish the crossover phase transition with a rapid 



Hydrodynamic afterburner for the CGC at RHIC 



4 



change of the thermodynamic variables from the first order phase transition. 

The statistical model and the blast wave model tell us that chemical freezeout 
temperature T ch is larger than thermal freezeout temperature T ch at the RHIC energies 
(see, e.g., discussions in Ref. [T7j). Chemical freezeout temperature obtained from the 
statistical model analysis is very close to the (pseudo) critical temperature of the QCD 
phase transition. This means that the hadron phase is almost chemically frozen due to 
strong expansion. For chemically frozen hadronic gas, one introduces chemical potential 
for each hadron [IB|. This modifies the chemical composition in the hadron phase 
compared with the chemical equilibrium hadronic gas and, consequently, changes the 
space-time evolution of temperature field "T"""| . 

1.2.2. Initial Condition Once initial conditions are assigned, one can numerically 
simulate the space-time evolution of thermalised matter which is governed by 
hydrodynamic equations. Usually, initial conditions are parametrised based on some 
physical assumptions. Transverse profile of the energy /entropy density is assumed to 
be proportional to the local number of participants p pa rt, the local number of binary 
collisions p co n or linear combination of them. Initial transverse flow is usually taken 
to be vanishing. In either full three dimensional (3D) hydrodynamic simulations 
or hydrodynamic simulations of longitudinal expansion with cylindrically symmetric 
geometry, one needs to parametrise also longitudinal profile of initial conditions 
for energy density, baryon density and four fluid velocity. The initial longitudinal 
shape is chosen so as to reproduce the final (pseudo)rapidity distribution of hadrons. 
Unfortunately, the shape is not uniquely determined. Two completely different initial 
conditions can end up almost similar rapidity distribution. See, e.g., Ref. [20J - 

On the other hand, one can introduce model calculations to obtain the initial 
condition of hydrodynamic simulations. Event generators can be used to obtain the 
energy density distribution at the initial time. Recently, the SPheRIO group employs 
an event generator NeXus and takes an initial condition from this model in the event- 
by-event basis [21]. The resultant energy density distribution in the transverse plane has 
highly bumpy structures J2T] 122]. Smooth initial conditions used in the conventional 
hydrodynamic simulations are no longer expected in one event. Another important 
example which is relevant at very high collision energies is an initial condition taken 
from the Colour Glass Condensate picture [23]. This will be discussed in details in 



1.2.3. Freezeout Conventional prescription to obtain the invariant momentum spectra 
from the hydrodynamic simulations is to employ the so-called Cooper-Frye formula [24*] 



where = (E, p) and d{ are, respectively, a four momentum and the degree of freedom 
of a particle under consideration. Integral is performed on the freezeout hypersurface 
E. Here da is an element of freezeout hypersurface on E. One usually assumes 



Sec. 2. 




d 3 p (2tt) 3 h{x) exp[(p"u v (x) - m(x))/T{x)} ± 1 



(3) 
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£ = S(T = T th or e = e th ), where T th (e th ) is thermal freezeout temperature (energy 
density). T th or e th is an adjustable parameter for reproduction of the slope of 
spectrum. Information on the hydrodynamic simulations is taken through da, w M , T 
and /ij at a space-time point x. Chemical potential for each hadron ^ is not vanishing 
in general even in vanishing baryon density due to T ch > T th . Although the chemical 
potential is important to reproduce the particle ratios as well as the shape of the particle 
spectra, this is often neglected for simplicity in the literature. 

The physical picture described by the Cooper-Frye formula is sometimes called 
"sudden freezeout" since the mean free path A is suddenly changed from zero to infinity 
through a thin freezeout hypersurface. Instead of using this, one can use a hadronic 
cascade model to describe the space-time evolution of hadrons j^o"] 125] . The mean 
free path among hadrons is finite and depends on hadronic species. Hence, one can 
describe a continuous freezeout picture through hadronic transport models. Note that 
continuous particle emission can be considered within the hydrodynamics [21] • Adoption 
of hadronic transport models after hydrodynamic evolution of the QGP liquid could 
refine the dynamical modeling of relativistic heavy ion collisions. However, it is not so 
easy to connect them in a systematic and proper way. One has to use Eq. (j^J again for 
the boundary between the QGP liquid and the hadron gas. However, Eq. (j3J) describes 
not only out-going particles p^da^ > but also in-coming particles p^da^ < 0. If one 
uses hydrodynamic simulations as inputs for hadronic cascade models, one must discard 
the in-coming particles, which violates the energy-momentum conservation. In order 
to remedy this problem, one need to solve hydrodynamic equations together with the 
appropriate kinetic equations [2%] . 

1 . 3. Hydrodynamic Results for v 2 

We show comparison of the hydrodynamic results with the experimental data for V2(pt), 
^2(77) at the RHIC energy and the excitation function of i^. 

1.3.1. V2{pr) With a help of assuming the Bjorken flow [2H] for the longitudinal 
direction, one can solve the hydrodynamic equations only in the transverse plane at 
midrapidity. Systematic studies based on this (2+l)-dimensional hydrodynamic model 
are performed in Ref. ]2J. For the EoS, complete chemical equilibrium is assumed 
for both the QGP phase and the hadron phase, dependences of t>2 for pions and 
protons from this model [H] are compared with the STAR data [30] in Fig. |2j (left). By 
employing the EoS with phase transition, the hydrodynamic model correctly reproduces 
v 2 (pt) and its mass-splitting behavior below px — 1 GeV/c. On the other hand, 
V2(pr) for (anti)protons from the resonance gas model does not agree with the data. 
Although the reason for the difference of the result between these two EoS models is 
not so clear, the experimental data favors the QGP EoS. Due to the assumption of 
chemical equilibrium in the hadron phase, this model does not reproduce particle ratio 
and spectra simultaneously. It is of importance to study whether the agreement with 
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Figure 2. (left) Transverse momentum dependence of v-i for pions and (anti)protons 
|Hj are compared with the STAR data [3D. in minimum bias Au+Au collisions at 
t/Fnn = 130 GeV. EoS Q stands for the equation of state with first order phase 
transition between the massless free parton gas and the resonance gas. On the other 
hand, EoS H stands for the equation of state for the resonance gas. (right) V2 for 
charged hadrons as a function of pseudorapidity 77 in Au+Au collisions at ^snn = 130 
GeV m321- PCE means the EoS of partial chemical equilibrium, whereas CE means 
the EoS of chemical equilibrium. Figure is taken from Ref. |19| . 

the experimental data still holds even when the assumption of chemical equilibrium in 
the hadron phase is abandoned j!9j . 

1.3.2. v 2 (r)) One needs a full 3D hydrodynamic simulation to obtain the rapidity 
dependence of v 2 since either cylindrical symmetry or Bjorken flow |2H| keeps us 
from obtaining this quantity. First analysis of v 2 (j]) at RHIC based on the full 
3D hydrodynamic model is performed in Ref. |31j . After that, the effect of early 
chemical freezeout on the elliptic flow is also investigated Figure El (right) 
shows V2{rj) for charged hadrons in Au+Au collisions at ^/Sjvv = 130 GeV [TJ 132] . 
Here the initial condition for the energy density is so chosen as to reproduce the 
pseudorapidity distribution of charged hadrons. Hydrodynamic results are consistent 
with the experimental data only near the midrapidity while hydrodynamics overpredicts 
the data in the forward/backward rapidity regions. Multiplicity is not so large in the 
forward/backward rapidity regions, so equilibration of the system tends to be spoiled. 

1.3.3. Excitation function Figure El shows the excitation function of v 2 compiled by 
the STAR Collaboration [23] ■ Hydrodynamic results presented in this figure are based 
on the same model discussed in Sec. 11.3.11 Data points continuously increase with the 
unit rapidity density at the AGS, SPS and RHIC energies. However, the hydrodynamic 
response is almost flat or slightly decreases with the multiplicity. The data points seem 
to reach the "hydrodynamic limit" for the first time at the RHIC energy. It will be 
interesting to see the same plots including the upcoming LHC results: If the system 
produced at LHC obeys the hydrodynamic picture, this experimental "curve" will bend 
at (l/S)dN/dy ~ 25-30. 
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Figure 3. Excitation function of vije. Figure taken from Ref. .''>.'? . 



The deviation between the hydrodynamic results and the data plots below 
(l/S)dN/dy = 25 reminds us the pseudorapidity dependence of elliptic flow in Fig. EI 
(right). The deviation might come from a common origin [33]: the small multiplicity 
both in forward rapidity region at the RHIC energy and at midrapidity at the lower 
collision energies could cause the partial thermalisation. In the low density limit, v 2 is 
actually proportional to the number density [35] ■ The shape of v 2 (r]) data in forward 
rapidity region looks similar to that of the pseudorapidity distributions jSH]- Similarly, 
data plots of excitation function increase almost linearly with the particle density. These 
results suggest that thermalisation is not achieved completely in forward rapidity region 
at the RHIC energy and at midrapidity at the SPS energies. 

Figure El shows the excitation function based on the hydro+cascade model [2E] with 
Bjorken longitudinal flow [21]. In the hadron phase, a hadronic event generator RQMD 
is employed in these calculations. Contrary to the conventional hydrodynamic models, 
freezeout processes are automatically described by this model without any adjustable 
parameters. The excitation function in the case of the latent heat ~ 0.8 GeV/fm 3 
linearly increases with the multiplicity, which is consistent with the experimental data. 
It should be noted that v 2 (pt), its mass dependences and particle spectra/ratio at 
midrapidity are also reproduced by this hybrid model [2TP . 

1.4- Summary and Outlook 

The most successful approach based on the hydrodynamic model to elliptic flow in 
relativistic heavy ion collisions is the hybrid one in which the QGP phase is described 
by the ideal hydrodynamics while the hadron phase is described by a hadronic cascade. 
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Figure 4. Excitation function from the hydro+cascade model. Figure taken from 
Ref. 

The Bjorken scaling solution [2H] is assumed in the current hydro+cascade model. This 
means that current hybrid models are available only near midrapidity. Therefore it is 
desired to develop a model in which a full 3D hydrodynamic simulation is combined 
with a hadronic cascade model. From agreement of excitation function between the 
hybrid model and the experimental data at midrapidity, the 3D hybrid model is 
expected to reproduce the pseudorapidity dependence of elliptic flow. It should be 
emphasised again that the hybrid model has its own problem on the violation of energy 
momentum conservations between the QGP liquid and the hadron gas. A simulation 
which incorporates a proper treatment at the boundary between the QGP phase and 
the hadron phase is now an open problem. 

2. Hydrodynamic Model with a CGC Initial Condition 

2.1. Toward a Unified Description of Relativistic Heavy Ion Collisions 

As already discussed in the previous section, one of the important findings at RHIC is 
that the hydrodynamic approach to the description of elliptic flow works remarkably well 
when early thermalisation time and the QGP EoS are assumed |H] . Hydrodynamics 
should work only in the situation that the mean free path A = 1/erp is much smaller 
than the typical size of the system L. This suggests that the particle density p becomes 
very large at RHIC. 

In addition to this, there are many important findings in heavy ion collisions at the 
RHIC energy. One of them is jet quenching |3Z|, i.e., suppression of high pt hadrons in 
the single particle spectra and disappearance of away-side peak in the di-hadron spectra. 
These phenomena are not found in dAu collisions [SHI EH EDJ \H\ : the yield of high pt 
hadrons roughly scales with the number of binary collisions N co \\ and the away-side peak 
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appears in the di-hadron spectra. So one can conclude that high p T partons produced 
in the initial hard scattering interact with the highly dense medium in the final state 
and that these partons lose their energies during traversing medium. 

So the current RHIC data strongly suggest that the parton density created in heavy 
ion collisions at RHIC is dense enough to cause both large elliptic flow of bulk matter and 
large suppression of high-p^ hadrons. What is an origin of this dense matter at RHIC? 
The bulk particle production in high energy hadronic/nuclear collisions is dominated 
by the small x modes in the nuclear wave function, where x is a momentum fraction 
of the incident particles. Hence, an origin of the large density could be traced to the 
initial parton density at small x inside the ultra-relativistic nuclei before collisions. It 
is well known that the gluon density increases rapidly with decreasing x until gluons 
begin to overlap in the phase space where nonlinear interaction becomes important [42J. 
These gluons eventually form another type of extremely dense matter, the Colour Glass 
Condensate (CGC) jl^j. This phenomenon is characterised by a "saturation scale" Q 2 S 
which is identified with the gluon rapidity density per unit transverse area. When Q s 
is sufficiently larger than the typical hadronic scale (Q s 3> Aq C d), the strong coupling 
constant becomes weak (a s (Q s ) < 1). Moreover, the gluon occupation number in the 
wave function becomes huge, ~ l/a s (Q s ) 3> 1. Therefore the CGC can be studied by an 
weak coupling classical method also known as McLerran-Venugopalan (MV) model jUj. 
The calculations based on the CGC picture have been compared to various RHIC data. 
Remarkably, the CGC results on the global observables, e.g., the centrality, rapidity and 
energy dependences of charged hadron multiplicities agree with the RHIC data under 
the assumption of parton-hadron duality 

In relativistic heavy ion collisions, the CGC could be a seed of the QGP 
So hydrodynamic initial conditions can be calculated from the two CGC collisions. 
One needs a non-equilibrium model which describes equilibration from the produced 
gluon distribution to a local thermal distribution. This is, however, beyond the present 
work. In this work, the rapidity, transverse position and centrality dependences of the 
number of produced gluons are identified with the hydrodynamic initial condition of a 
thermalised QGP. 

Some of the problems which are inherent in a particular approach can be removed. 
For instance, one conventionally parametrises initial conditions in the hydrodynamic 
simulations shortly after a collision of two nuclei by hand. Therefore it is desired to take 
an initial condition which is obtained by a reliable theory. On the other hand, most of 
the calculations based on the CGC do not include final state interactions. According to 
estimation of initial parton productions from a classical lattice Yang-Mills simulation 
based on the CGC jSlllHj, the transverse energy per particle is roughly E T /N ~ 1.5 
GeV, while the experimental data yields Et/N ~ 0.6 GeV jlHl- Moreover, elliptic 
flow from a classical Yang-Mills simulation is inconsistent with RHIC data [SD]- Those 
discrepancies, which are mainly due to the lack of collective effects, will be removed by 
introducing further time evolution through, e.g., hydrodynamics. For the calculation of 
parton energy loss, one needs time evolution of parton density. Bjorken expansion [2H| 
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is assumed for the evolution of parton density in almost all work except Refs. |oT] 152] . 

We have already developed a unified model (CGC+Hydro+Jet model) in which 
hydrodynamic evolution, CGC initial conditions, and the production and propagation 
of high p T partons are combined [23] • Only the soft sector, i.e., the production of gluons 
in CGC collisions and the hydrodynamic afterburner for the CGC initial condition 
is discussed in this paper. For further details of the model and the results for high 
Pt regions such as centrality dependences of the nuclear modification factors and the 
azimuthal correlation functions, see Ref. [23] . 

2.2. Gluon Production from a CGC picture 

We employ the kx factorisation formula for the production of gluons from CGC collisions 
along the line of work by Kharzeev, Levin and Nardi (KLN) |UJ. The number of 
produced gluons in the /^-factorisation formula is given by [121 E3 EH ES] 

dN g 4ir 2 N c rd 2 p T 
d 2 x±dy N 2 — 1 J p\ 

x 4> A {xi, x_ L )0 jB (x 2 , (pt ~ k T ) 2 ] x±), (4) 

where Xi j2 = px exp(±y) / y/s , and y and px are a rapidity and a transverse momentum 
of a produced gluon. Running coupling constant a s is evaluated at the scale Q 2 = 
m&x(k T , (p T — k T ) 2 ). The unintegrated gluon distribution (p is related to the gluon 
density of a nucleus by 

rQ 2 

xGa(x,Q 2 ) = / d 2 kTd 2 x±(pA(x, /c^;x_|_). (5) 

Motivated by the KLN approach, we use a simplified assumption about the unintegrated 
gluon distribution function: 

^,4;xJ =(^#' k , T l Q " (6) 

where Cp = (N 2 — 1)/(2N C ). We introduce a small regulator A = 0.2 GeV/c in order 
to have a smooth distribution in the forward rapidity region \y\ > 4.5 at RHIC. Other 
regions are not affected by introducing a small regulator. The above distribution depends 
on the transverse coordinate x^ through Q 2 . Saturation momentum of a nucleus A in 
A + B collisions is obtained by solving the following implicit equation with respect to 
Q s at fixed x and x^ 

2tt 2 

Q 2 s (x,x ± ) = — a s (Q 2 s )xG(x, Q 2 s )p^ vt (x±) } (7) 

L> F 



where 



Pirt(xx) = T A ( X± ) {l - [1 - a^ N T B ( X± )/B] B } . (8) 
The number of participants is given by 

iVpart = / d 2 X ± ( P ^ rt (x ± ) + pf art (x ± )). (9) 
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We take a simple perturbative form for the gluon distribution in a nucleon 

xG(x,Q 2 ) = K\n ( Q2 J~ A ) x-\l - x) n (10) 

V A QCD / 

where A = Aq CD = 0.2 GeV. K is used to control saturation scale in Eq. (J7|) [56j. We 
choose K = 0.7 for A = 0.2 so that the average saturation momentum in the transverse 
plane yields {Q 2 {x = 0.01)) ~ 2.0 GeV 2 /c 2 in Au+Au collisions at impact parameter 
b = at RHIC. Similar to the KLN approach, x~ x dependence of the saturation scale is 
motivated by the Golec-Biernat-Wiisthoff model ^7\. The factor (1 — x) n shows that 
gluon density becomes small at x — > 1. n usually depends on the scale Q 2 . Here we 
take n = 4 as in the KLN approach [35] . 

We obtain the rapidity distribution for produced gluons at each transverse point xi 
by performing the integral of Eq. (0J) numerically. The transverse energy distribution 
dEx/dy is also obtained by weighting the transverse momentum of gluons i n the 
integration with respect to in Eq. (jH). We cut off the integral range of p? in Eq. (JU) 
since only the low p? partons are assumed to reach the local thermal equilibrium. We set 
the cut-off momentum as pr.cut — 3 GeV/c which corresponds roughly to the maximum 
saturation scale at x = 0.01 at the origin = in central Au+Au collisions at RHIC. 

In order to obtain initial conditions for hydrodynamics, one needs a non-equilibrium 
description for the collisions of heavy nuclei. Leaving this problem for the future 
work, we assume that the system of gluons initially produced from the CGC reaches a 
kinematically and chemically equilibrated state at a short time scale. 

There are two ways to provide initial conditions from the CGC, i.e., matching 
of number density (IC-n) and matching of energy density (IC-e). Here, we employ 
the former prescription, i.e., IC-n. For the result from IC-e, see Ref. [23J. Assuming 
Bjorken's ansatz y = r] s [2H| where rj s is the space-time rapidity r] s = (1/2) ln(x + /x~), 
we obtain the number density for gluons at a space-time point (r , xj_, r] s ) = (r , x) from 
Eq. © 

n ^ x) = 7^x- ± - 

Number density n for the massless free parton system can be written: 

n= {^d q + d g )^T\ (12) 

where d g = 2N c N s N f = 36, d g = 2{N 2 - 1) = 16 and C(3) = 1.20206. Thus, we obtain 
the initial temperature field in the three-dimensional (3D) space from the number density 
(IC-n): 

T(r °' x) = i^3C(3r| • (13) 

When the temperature at (r , x) is below the critical temperature T c = 170 MeV, the 
all thermodynamic variables in the cell are set to zero. 
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Figure 5. Pseudorapidity distributions of charged hadrons in Au + Au collisions at 
y/sWN= (a) 130 and (b) 200 GeV are compared to the PHOBOS data [SHj- Figure 
taken from Ref. 
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Figure 6. Mean transverse momenta for pions, kaons and protons as a function 
of -/V part |59| . Here contribution only from hydrodynamic components is taken into 
account. T th = 100 MeV is used for all centralities. Figure taken from Ref. 



2.3. Results 

In Fig. pseudorapidity distributions of charged hadrons in Au + Au collisions at 
both ^snn = 130 and 200 GeV are compared with the PHOBOS data [HE]- Impact 
parameters in each panel are, from top to bottom, 2.4, 4.5, 6.3, 7.9, and 9.1 fm (2.5, 4.4, 

6.4, 7.9, and 9.1 fm) for ^/s^w = 200 (130) GeV. These impact parameters are evaluated 
from the average number of participants at each centrality estimated by PHOBOS [58j. 
It should be also emphasized that it is not easy to parametrize such initial conditions 
which reproduce the data with the same quality as the CGC initial conditions presented 
here. 

Mean transverse momenta (pt) for pions, kaons and protons as a function of A^ part 
are compared with the PHENIX data [SH] in Fig. |H1 Although our results are slightly 
smaller than the data in central and semicentral regions, the overall trend is consistent 
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with the data. For pions, the semihard spectrum starts to be comparable with the 
soft spectrum around px = 1.5-2.0 GeV/c. We reproduce the px spectra for pions by 
including contribution from quenched jets. So the deviation for pions can be filled by 
the semihard spectrum. While semihard components for kaons and protons are very 
small in low and intermediate pt regions. A little more radial flow is needed to gain the 
mean transverse momentum in central collisions. 
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Figure 7. (Pseudo)rapidity distribution of (a) gluons produced from CGC collisions, 
(b) partons at initial time, and (c) final total hadrons in Au + Au collisions at */sjvjv = 
200 GeV. 

Both the KLN approach and the hydrodynamic model with a CGC initial condition 
reproduce the centrality and rapidity dependences of particle distributions. Then, 
what is the role of hydrodynamic evolution in comparison with the KLN approach? 
The main difference between these two approaches is whether final state interactions 
are taken into account. The initial transverse energy per particle is estimated to be 
Er/Ng ~ 1.6 GeV at y = 0. It becomes Et/N 9 = 1.0 GeV after assuming a thermal 
state. The effect of the hydrodynamic afterburner is to reduce the transverse energy 
per particle due to pdV work. We find that (dEx / dy) / (dN / dy)\ y=0 = 0.54 GeV and 
(dE T /dr])/(dN/dr])\ v=0 = 0.56 GeV. 

2-4- Summary and Outlook 

Our approach is a first step toward a unified dynamical modeling of relativistic heavy 
ion collisions. 

The main goal has been developing a consistent dynamical framework for the space- 
time evolution of both bulk matter and hard jets. In order to accomplish this, much 
remains are to be done. 

• More sophisticated wave function should be used to calculate the production of 
initial gluons. At midrapidity, results from the classical Yang-Mills equation can 
be used. In forward rapidity region at the RHIC energy or even at midrapidity at 
the LHC energy, one needs a quantum evolution. 

• As discussed in the Sec. 1, the hadron phase should be described by a hadronic 
transport model. This will improve v% in forward rapidity regions in the 
hydrodynamic result. 
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• EoS from the recent lattice QCD simulations should be used in the hydrodynamic 
simulations. 

• Chemical non-equilibrium process will be taken into account by solving the rate 
equation for the number of quarks and gluons together with the hydrodynamic 
equations [HI]]- The amount of energy loss of partons depends on a species of both 
probe parton and medium parton. So the correct chemical composition in the 
QGP phase will lead to reproduction of nuclear modification factors without any 
free parameters. 
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